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ABSTRACT 

Among various approaches for indirect detection of dark matter, synchrotron emission due to secondary 
electrons/positrons produced in galactic WIMPs annihilation is raising an increasing interest. In this paper 
we propose a new method to derive bounds in the m^-{af^v) plane by using radio continuum observations 
of Messier 33, paying particular attention to a low emitting Radio Cavity. The comparison of the expected 
radio emission due to the galactic dark matter distribution with the observed one provides bounds which are 
comparable to those obtained from a similar analysis of the Milky Way. Remarkably, the present results are 
simply based on archival data and thus largely improvable by means of specifically tailored observations. The 
potentiality of the method compared with more standard searches is discussed by considering the optimistic 
situation of a vanishing flux (within the experimental sensitivity) measured inside the cavity by a high resolution 
radio telescope like ALMA. Under the best conditions our technique is able to produce bounds which are 
comparable to the ones expected after five years of Fermi LAT data taking for an hadronic annihilation channel. 
Furthermore, it allows to test the hypothesis that space telescopes like Pamela and Fermi LAT are actually 
observing electrons and positrons due to galactic dark matter annihilation into leptons. 
Subject headings: 



1. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs) still repre- 
sent the most natural candidate for Dark Matter (DM) whose 
detection is one of the main purposes of astrophysical and 
cosmological observations nowadays (iBertone et al. 20051 
iKomatsuet al. 2008). 

From a theoretical point of view, models of WIMPs arise 
straightforwardly in supersymmetric extensions of the Elec- 
troweak Standard Model as the lightest supersymmetric par- 
ticle. These candidates can annihilate in couples to produce 
ordinary particles as final states which, in principle, can be 
detected. 

Concerning the produced electromagnetic radiation, 
gamma-rays provide the most promising opportunity due 
to the very low attenuation in the interstellar medium, 
and to its high d et ection efficiency (see for example 
lJungmanetal. 19951 iBergstrom 20001 IBertone et al. 20051 
for a review of this extensively studied issue). In this respect 
the activity of the recently launched Fermi Large Area 
Telescope provides an exciting chance to have a large amount 
of data on the gamma-ray sky between 20 MeV and 300 GeV, 
inclu ding the largely u nexplored energy window above 10 
GeV (lAbdo et al. 20081) . 

Nevertheless, high energy electrons and positrons, 
once produced in the annihilation cascade of WIMP 
pairs, interact wi th gas, radiation an d the magnetic field 
in the galaxv (iBaltz & Edsio 19981 iDonato et al. 2003L 
iHooper&Silk 2004 . During their propagation in the 
galactic medium they release secondary radiation in th e radio 
and X-ray bands (Blasi et al. 2002, T asitsiomi et al. 20031 
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l Aloisio et al."200l iBaltz & Wai 2004 IColafrancesco 2005 



IColafrancesco 200a, IBe rgstrom et al. 2006 

Keds & Ullio 2008, Heltema & Profumo 2008, 

IZhang & Si gl 2008), giving rise to a further signal of 
the DM annihilation that i s, in principle, detectable. 

In lBorrielloetal. 20091 some of the authors of the present 
paper have developed a semianalytical approach to com- 
pute the synchrotron emission from the Milky Way (MW) 
WIMP pairs annihilation. The comparison of the the- 
oretical DM radio signal from our Galaxy with the ob- 
served emission, (where only the CMB is modelled and re- 
moved) allows to draw an exclusion plot in the ■m^-{aAv) 
plane. For this purpose one can use the code described in 
Ide Oliveira Costa et al. 200* 81 where most of the whole sky 
radio observations in the range 10 MHz- 100 GHz are col- 
lected and a scheme is described to derive interpolated, CMB 
cleaned sky maps at any frequency in this range. This result, 
however, is strongly affected by two problems: 

i) The difficulties in removing the foreground contribu- 
tion from the radio continuum. Such a removal would 
allow, at least, an important refinement of the m^- 
{(Jav) bounds or, even better, the detection of DM. Let 
us recall that different approaches to the foreground re- 
moval at the WMAP frequencies produced very differ- 
ent results, exhibiting a residual signal (which has been 
interpreted as a signatur e of the DM anni hilation) as 
in IFinkbeiner 2004 and Hoo per at al. 20071 or no addi- 
tional feature, as in a mor e recent analysis performed 
by the WMAP team itself (iGold et al. 2008h . 

ii) The fact that the measured quantities are averaged by 
the integration along the line of sight which prevents 
a proper estimate of the local values of both galactic 
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Fig. 1 . — Radio continuum emission from M33 at A = 3.6 cm (left) and A = 20 cm (riglit), witli apparent B-vectors of polarized intensity superimposed. The 
beam areas are shown in the left corner of each figure. Contour lev els are 0.5, 1, 2, 3, 4, 8, 16, 32 and 64 mJy/beam for the 3.6 cm and 0.3, 0.6, 0.9, 1.2, 1.8, 2.4, 
3.6, 4.8, 9.6 and 19.2 mJy/beam at 20 cm ITabatabaei et al. 2007bl) . 



magnetic field and radio continuum emission. The high 
level of uncertainty on the local values of the MW mag- 
netic field propagates on the results after the integration 
along the line of sight (an order ten kpc long path) is 
performed. 

To overcome problems ii) we will reapply such an anal- 
ysis to the nearby, close to face-on, Scd galaxy Messier 
33 (M33), taking advantage of the thinness of the re- 
gion source of the signal. This galaxy is the third 
largest in the Local Group and is located at a dis- 
tance of 840 kpc from us dFreedman et al. 19911) . It 
has been extensively studied at all wavelengths: radio 
continuum (Tabatabaei et al. 2007b, Tabatabaei et al. 2007c, 
ITabatabaei et al. 2008), far infrared (Tabatabaei et al. 2007a), 
near IR fCioni et a l. 2007). optical (.Mass _ev_elal._^006J, 
UV dC iani et al. 1981 iGordonetal. 1999I) . " and X-ray 
jHaberl & Pietsch 200l1). Its kinemati cs has also been the 
subject of many studies dCarignan 20061) . M33 has an inclina- 
tion of ~ 56° and due to its rather flocculent nature, it presents 
well defined arm and interarm regions. Some of the latter re- 
gions can be characterized by a particularly low level of ra- 
dio emissivity (hereafter denoted for simplicity Radio Cavi- 
ties), even though placed not far from the galactic center, and 
thus where possibly a non negligible magnetic field is at work. 
Hence, the poor radio emission would be simply due to a lack 
of astrophysical electrons/positrons. In this framework, the 
presence of a halo of WIMPs would inject in any case a flux 
of electron/positrons. These particles, while diffusing in the 
galaxy and propagating in the magnetic field, have to produce 
a radio signal compatible with the observed one. The require- 
ment that the DM induced radio signal coming from a particu- 
lar Radio Cavity of M33 does not overcome the observed one 
straightforwardly fixes bounds in the m^-{(jj^v) plane. 

Analyses devoted to the study of radio emission by exter- 
nal galaxies are alread y pr esent in literature (se e for example 
iTasitsiomi et al. 20031 and IColafrancesco 200 ^, even if not 
focussing on Radio Cavities. It is worth noticing that the main 



advantage of considering a small region rather than the whole 
galaxy is that this approach allows to keep the different astro- 
physical hypothesis much more under control. For the pur- 
poses of the present work, M33 provides an ideal case since 
it is near, it has a small bulge-to-disk ratio, and therefore, it 
is possible to trace the inter-arm regions also relatively close 
to the central and densest regions. Moreover, due to its rela- 
tively small inclination (e.g. compared to M31), it allows to 
minimize the effects of the integration of the signal along the 
line of sight. 

Unfortunately we cannot solve the first problem. The ra- 
dio emission of galaxies like M33 has three components: the 
clumpy thin disk (~ 200 pc scale height) with thermal and 
synchrotron emission from star-forming regions and super- 
nova remnants, the smooth thick disk (^ 1 kpc) with mostly 
synchrotron emission by young electrons originating from the 
thin disk, and the smooth halo (unknown extent) with syn- 
chrotron emission by old electrons from the thin disk. The 
relative contributions may vary from galaxy to galaxy, but in 
general the radio from the thic k disk dominates, as visible in 
edge-on galaxies dBeck 20091) . The emission from the dark 
matter halo will mix with the thick disk emission but the mix- 
ing ratio is of course unknown. For this reason we preferred a 
conservative approach not subtracting any possible astrophys- 
ical component. 

The paper is organized as follows. In Section|2]we describe 
the main characteristics of M33 Radio Cavity, whereas in Sec- 
tion 13 we compute the expected radio signal coming from it. 
Finally, in Sec. |4]we derive our constraints for WIMPs in the 
m^-{<TAv) plane and give our conclusions. 

2. THE RAD/0 CAV/rF IN M33 

In the present analysis, we use the observations of 
M33 at A = 3.6 cm and A = 20 cm described in 
ITabatabaei et al. 2007bl and reported in Fig. [l] 

A detailed description of the data acquisition reduction and 
analysis can be found in the above quoted reference. Here we 
shall just point out a few relevant facts. The data at A = 3.6 
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Fig. 2. — Emission from M33 at different wavelengths. From the upper left 
clockwise: radio continuum at 3.6 cm, radio continuum at 20 cm, far infrared 
at 8.0 fim (data from Spitzer telescope), and Near UV at 2300 A. The Radio 
Cavity is marked by a circle. 



cm were obtained using the 100-m Effelsberg radiotelescope 
and the final map has a HPBW of 80" and an r.m.s noise of 
~ 220 /iJy/beam. On the other hand the data at A = 20 cm 
were obtained using the VLA-D array and corrected for lack 
of flux, due to missing short spacing, using single dish 20 cm 
observations obtained at the Effelsberg radiotelescope. The 
resulting map has an HPBW of 51" and an rm.s noise of ~ 
70 ^Jy/beam. 

For the purpose of this work, the most noticeable feature, 
marked as a circle in Fig. |2] is the interarm cavity located at 
RA ~ lh33'"24.0'' and Dec ~ 30° 35' 39.0", i.e. at ~ 0.14° 
SE of the bright nucleus (2. 1 kpc from the center). This cavity 
has an area of '--^ 5 x 10'^ arcsec^, and is therefore well resolved 
at all wavelengths. 

Despite the small optical and radio flux coming from the 
cavity, there is no reason to believe the magnetic field there 
to be significantly smaller than the one in the surrounding re- 
gions. Therefore, the cavity is the ideal place to look for a 
DM signal, since there the ratio between the DM radio signal 
and the astrophysical foreground reaches its maximum value. 

3. DM SYNCHROTRON SIGNAL 

Let us first assume that M33 is perfectly face-on (later we 
will correct for the angle it forms with the line of sight) and 
consider a cartesian galactocentric coordinate system (x, y, z) 
with the z-axis pointing in our direction. If is the syn- 
chrotron emissivity due to DM annihilation, the radio flux 
coming from an elementary volume dV within M33 is given 
by dl^ — j^dV/AnD'^, where D is the distance of M33 from 
us. 

Thanks to the face-on approximation and the large value 
of D, we can simply write dV = D^dVldz and consider the 
integration along the line of sight as a simple integration along 
the z-axis. At last, the thinness of the integration volume 
allows us to take into account the angle 6 the galaxy forms 
with the line of sight by simply multiplying the flux by cos 6. 
The result is 



di^ cos e 



where a — (x/D) cos 9 and (3 ^ y/D are two visual angles 
and where we chose the y-axis to be coincident with the M33 
major axis. 

We perform the integration of Eq. [T] over the range z G 
[—2 kpc, 2 kpc]. This choice does not sensibly affect the re- 
sults due to the exponential decrease of the DM synchrotron 
emissivity away from the galactic disk. 

An analytical e xpression of can be found in 
iBorriello et al. 20091 in terms of the DM density profile, the 
magnetic field and the interstellar radiation field (ISRF). The 
remaining parts of this section are therefore devoted to the 
parametrization of the Radio Cavity in M33. 

3.1. Dark Matter density profile 

It is not yet possible to uniquely determine the DM den- 
sity of M33, nonetheless both a cored and a spiked profile 
ca n be deduced fi tting the galaxy rotation curve (see Fig. 5 
of Corbelli 200 3il). Th e cored case is well represented by a 
Burkert (Burker t 1995h density profile 
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while the spiked profile is well descr ibed by and a NEW den- 
sity distribution dNavarro et al. 1996h 
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We consider two of the various fits obtained in lCorbelli 2003L 
Their defining parameters are: 
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Profi les with a centr al slope steeper tha n r ^ like the Moore 
profile (iMooreetal. 19 99) are excluded dCorbelh 2003) . 

3.2. Magnetic Field 

Following iTabatabaei et al. 20081 the equipartition total 
magnetic field strength in the radio cavity is estimated as 
7. 1 ± 0.5 fjG. An exponential decrease is assumed to describe 
the field along the z axis: 



B = Boe- 



(2) 



jy(x,y,z)dz , 



(1) 



where zq is the scale height of the magnetic field. Observa- 
tions of edge-on galaxies give a sc ale height o f about 1.8 kpc 
for the synchrotron emission (see iBeck 20091) . This leads to 
a magnetic field scale height of zq ~ 7 kpc, assuming a non- 
thermal spectral index of ^ 1 (see e.g. Klein et al. 1982). 

In the case in which the equipartition hypothesis were not 
valid the situation would be more involved and two cases can 
be taken into account: 

(a) If the field strength is constant across the cavity, the cos- 
mic ray electrons (CRE) have to be deficient to explain the 
low radio emission. B can hardly be constant with height, 
too, but it could decrease even slower than with a scaleheight 
of 7 kpc if CREs suffer from strong energy losses (which is 
probable). We estimated a 7.5 pG field strength in this case 
and a scaleheight that could vary from 5 to 10 kpc; 

(b) If however the CRE density is about constant across the 
cavity, the field there would be smaller than 7 pG and, if it is 
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Fig. 3. — Constraints in the m;^-((T^ti) plane from radio observation of the Radio Cavity of M33 at 3.6 and 20 cm. Full annihilation into bb couples is assumed. 
The upper panels refers to a NFW DM density profile, the lower ones to a Burkert profile. Three cases are considered: (1) the equipailition hypotesis is verified; 
(2a) it is not and B is constant across the cavity; (2b) it is not and CRE density is constant acros s the cavit y. Solid lines refer to M33, while dot-dashed lines 
show the bounds obtainable - under the same physical hypothesis - for the MW (as in cBorriello et al. 20091) ). 



also constant with height, the S-scaleheight would be 2 times 
smaller. We estimated B ^ 3.2 /iG and zq = 3.5 kpc. This 
is however quite unlikely because CREs loose energy with 
height. 



3.3. Interstellar Radiation Field 

Once the high electrons are injected, they lose energy due, 
mainly, to two physical process: the synchrotron radiation we 
are interested in, and inverse Compton scattering off of the 
ISRF photons. Therefore, we need to characterize the ISRF 
in the cavity to properly subtract the amount of energy lost in 
the scattering wi th the photo ns. 

According to iDeul 19891 we assume for the ISRF in the 
diskof M33 ul°^{R) = 5.32 e--«/2.io ^j^ere x = 0.539 
eVcm^"^ is the MW's ISRF at the Solar System position 
jWeingartner & Draine 200l . The previous expression is ob- 
tained by setting the ISRF equal to x to a galactocentric ra- 
dius of 3.5 kpc. To deduce the scale height of the exponential 
decrease of the ISRF far from the disk, we rescale the MW 
field in the same way. A galactocentric distance of 2.11 kpc 
in M33 corr esponds to 5.12 kpc in the MW. Using the model 
described in lPorter & Strong 20051 it is possible to deduce that 
the MW's ISRF scale height at a radius of 5.12 kpc is equal 
to 3.72 kpc which, when rescaled to M33, gives 1.53 kpc. 

Putting things together we parameterize the ISRF in the Ra- 



dio Cavity as: 

C/™d(2) = C/oe-l^l/^°, 
with Uq = 1.05 eV cm^^, and ho = 1.53 kpc. 

4. RESULTS AND CONCLUSIONS 

Following the approach described above one can compute 
the expected radio emission due to electrons/positrons pro- 
duced in the annihilation in pairs of WIMPs forming the halo 
of M33, and in particular the signal coming from the Radio 
Cavity which cannot overcome the observed signal. Hence, 
by using the observations of M33 at A — 3.6 cm and A = 20 
cm, we obtained the exclusion plots reported in Fig. |3] and 
|4] Three cases are plotted for each wavelength and DM den- 
sity profile: The equipartition hypothesis is verified (case 1: 
B = 7.1 /iG, Zq = 7.0 kpc); The equipartition hypothesis is 
not verified and the field strength is constant across the cav- 
ity (case 2a: B ~ 7.5 fiG, zo = 10 kpc) or the CRE density 
is constant (case 2b: B = 3.2 /iG, zq = 3.5 kpc). Let us first 
consider the results shown in Fig. [3] For both wavelengths the 
bounds are compared to those corresponding to the MW (dot- 
dashed lines), de rived in the same particle physics hypothesis 
jBorriello et al. 2 009) in which DM particles fully annihilate 
into bb couples. The results are quite dependent on the DM 
profile chosen, due to the position of the radio cavity (quite 
near the galactic center). In fact, assuming a Burkert density 
profile one gets less constraining results. In fact, at 2.1 kpc 
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Fig. 4. — As in Fig. [3]but for DM particles fully annihilating into ^^/i couples. This time the benchmark case is represented by the Pamela/Fermi LAT/Hess 
favored region (as in lMeade et al. 20091) corresponding to the elliptical spot centered at about 10^ GeV and 10"^'^ cm^s~^. 



from the galaxy center the Burkert profile is about 2/5 of the 
NFW one, therefore we obtain six times weaker bounds from 
the cored case. 

It is quite encouraging to observe that, at least for a favored 
case of an equipartition magnetic field and for a NFW profile, 
this calculation provides comparable or stringent bounds on 
DM than the ones obtained from our Galaxy. Moreover the 
M33 radio data used were not tailored for this aim, and thus 
could be largely improved in terms of resolution (see below). 

Let us consider the case of very heavy WIMP directly an- 
nihilating into charged leptons as recently proposed in or- 
der to explain the positron excess observed by Pamela satel- 
lite (lAdriani et al. 2009) and the e^-e" excess seen by Fermi 
LAT (lAbdo et al. 20 08). This unexpected flux is made of 
electrons with energy greater that ^ 100 GeV. The syn- 
chrotron emission due to electrons of this energy is peaked 
around 100 GHz, quite far from the values we are consider- 
ing here. Fig. |4] shows the bounds for the case in which the 
DM particle fully annihilate into /i^/i^ couples. In this case 
the DM contributes only for 1 % to the observed flux at 20 
cm, while the situation improves at 3.6 cm, where the DM 
accounts for a order 10 % of the emission. Therefore, higher 
frequency observations would be needed to adequately con- 
sider this case. 

Among ground-based telescopes, highest resolution (0.03") 
and sensitivity (0.060 mJy) observations at high-frequencies 
{v ~ 100 GHz) will be achieved by the Atacama Large Mil- 
limeter Array (ALMA). In addition Rotation Measurements 



Synthesis (Heal d et al. 2000h at as many wave wavelengths 
as possible will lead to a more precise magnetic field study, 
hopefully removing the degeneracy of possibilities shown in 
Fig. OandlH Such a high angular resolution could improve 
our results. In fact if one can identify a small region inside 
the cavity from which one detects a vanishing flux (within the 
experimental sensitivity) this naturally results in a significant 
improvement of our bounds. To be quantitative let us consider 
ALMA observing at a 110 GHz in a compact configuration 
(50 antennas, 3.7 arcsec beam, 16 GHz bandwith). By using 
the ALMA Sensitivit y Calculator one gets a 0.01 mK surface 
brightness sensitivity for one day of exposure time. 

In Fig. |5]we report the exclusion plots corresponding to the 
above null detection assumption for a NFW DM density pro- 
file. The left panel shows the results obtained for annihilation 
into hh while the right one refers to ii^n^ ■ For the hadronic 
case we compare our bounds, derived for one day and one 
month of observation, with the foreca st (at 5a) correspo nding 
to five years of Fermi LAT data taking lBaltz et al. 20081 Inter- 
estingly radio observations perform better for heavy WIMPs. 
In the right panel our exclusion regions are superimposed 
to the eflipse fav ored by Pamela/Fermi LAT/Hess data (see 
iMeade et al . 2009 and references therein). 

To summarize, we derived the expected secondary radi- 
ation due to synchrotron emission from high energy elec- 
trons/positrons produced in DM annihilations in M33. The 
comparison of the expected DM induced emission with ra- 
dio continuum observations, especially focussed on a particu- 
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Fig. 5. — Exclusion plots corresponding to a null detection (within the experimental sensitivity) radio flux. The left panel refers to DM particles annihilating into 
b b quarks derived for one day and one month of observation. The dot-dashed line represents the forecast (at Scr) for five years of Fermi LAT data taking obtained 
by means of Galprop for optimized parameters (see Baltz et al. 2008 for details). The right panel refers instead to annihilation into fi'^ fJ.^ . The exclusion regions 
are superimposed to the ellipse favored by Pamela/Fermi LAT/Hess data (see .Meade e t al. 2009). 



lar low emitting Radio Cavity, allows to obtain bounds in the 
ra^-{ap^v) plane. By using not tailored archival data we are 
able to put bounds comparable with ones obtainable by full 
sky observations of the radio emission from our Galaxy. This 
suggests the possibility that focused observations of small and 
well known astrophysical object could gain a great predic- 
tive power. To be more quantitative on the potentiality of the 
method we have considered an optimistic scenario were a high 
resolution radio telescope like ALMA does not detect any flux 
(within the experimental sensitivity) in a region inside the ra- 
dio cavity. In this case the bounds we can derive on WIMPs 
parameters are better than the corresponding Fermi LAT fore- 
cast in the case of heavy WIMPs decaying into hadrons. As- 



suming instead annihilation into muons we could be able to 
rule out the region favored by Pamela/Fermi LAT/Hess in a 
single day of observation. 
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